First-principles Calculation of Effective Onsite Coulomb Interactions of 3d Transition Metals: 
Constrained Local Density Functional Approach with Maximally Localized Wannier Function 
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We present a new ab initio method for calculating effective onsite Coulomb interactions of itinerant and 
strongly correlated electron systems. The method is based on constrained local density functional theory for- 
mulated in terms of maximally localized Wannier functions. This scheme can be implemented with any basis, 
and thus allows us to perform the constrained calculation with plane-wave-based electronic-structure codes. We 
apply the developed method to the evaluation of the onsite interaction of 2)d transition-metal series. The results 
are discussed using a heuristic formula for screened Coulomb interactions. 
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Properties of itinerant and strongly correlated electron sys- 
tems are widely discussed with phenomenological Hamiltoni- 
ans such as Anderson-impurity 1 1] and Hubbard models 
These models embody essential aspects of correlated elec- 
trons, and can be solved numerically or sometimes analyti- 
cally for special cases. However, parameters in the models 
are often determined empirically so that the employed model 
should reproduce experimental results of interest. 

Recently, there has been growing interest in constructing 
model Hamiltonians from first principles ||^. The principal 
motivation for such a study is a so-called long-standing "be- 
yond LDA" problem. One of recent major approaches toward 
correlated electrons from first principles is combining density- 
functional theory within local density approximation (LDA) 
with dynamical mean-field theory (DMFT) |4]. The main idea 
of this method is to map a first-principles Hamiltonian onto a 
lattice-fermion model to which the DMFT method is applica- 
ble, so developments of tractable but rational mapping tech- 
niques are highly needed. 

An important notice in constructing ab initio model Hamil- 
tonians lies in an evaluation of an onsite interaction param- 
eter known as the Hubbard U . Transfer parameters can be 
easily obtained by representing a one-body ab initio Hamil- 
tonian with a spatially localized (or atomic) orbital, while the 
U parameter evaluated with a one-center Coulomb integral of 
the atomic orbital gives a rather large value relative to that de- 
duced from experiments 1 5, 6] . This is because a screening ef- 
fect of surrounding valence electrons is completely neglected 
in the evaluation of the value. A representative treatment for 
calculating U including the screening effect is a constrained 
approach [7,, ,8]. It gives an optimally screened U, because, 
in that calculation, the U value is obtained from a response of 
the system to a change in a local charge density, thus incorpo- 
rating the valence-electron-screening effect statically into the 
U calculation. 

An important point in performing the constrained calcula- 
tion is a choice of basis functions to define the charge density 
of localized electrons. The use of an atom-centered localized 
basis set like a linear muffin-tin orbital (LMTO) |9] facilitates 
this definition, so constrained calculations to date have been 



performed with LMTO-based electronic-structure codes. The 
method has been widely applied to investigations of effective 
Hubbard J7's of various correlated materials 1 10,.l LI . while it 
is known that the calculated effective onsite Coulomb inter- 
action includes an eiTor 1 11] arising from the atomic sphere 
approximation employed in the conventional LMTO calcula- 
tions fl^ . 

In this Letter, we present a new method for calculating 
an effective U ; that is, a constrained local density functional 
approach based on "maximally localized" Wannier functions 
(WF's) fisll . The maximally localized WF, obtained by min- 
imizing a spatial spread of WF's in real space, has a practical 
advantage that it can be computed with any basis functions. 
This property enables us to implement the constrained scheme 
in the plane-wave-based electronic-structure codes. In addi- 
tion, the maximally localized WF is regarded as a reasonable 
basis for constructing model Hamiltonians; the property of 
the WF being well-localized spatially ensures that the Hamil- 
tonian matrix represented by this basis is sparse and short- 
ranged. In the present study, we apply the developed method 
to the systematic analysis for effective onsite Coulomb inter- 
actions of 3d transition metals (3- 

An electronic structure of a transition metal consists of lo- 
calized d electrons forming a narrow band and itinerant sp 
electrons associated with a wider band. There are n local- 
ized electrons per atom for the ground state. An effective 
Coulomb interaction Ues between two d electrons in an atom 
is defined as an energy cost in the electron transfer process, 
2d" d"^^ + d"+^, which is written in terms of total ener- 
gies as 



C/cff = E{n + 1) + E{7i - 1) - 2E{n), 



(1) 



where E{n) is the ground-state energy of the system, while 
E{n ± 1) coiTespond to total energies for cases where there 
are n ± 1 localized electrons in a specific atom. It should be 
noted here that the number of onsite d electrons is regarded as 
an adiabatic parameter for the electronic state of the system; 
i.e., for the total-energy calculation, the charge density of the 
system is relaxed in the constraint that the occupation number 
of localized d orbitals of the specific atom is kept at a given 
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value. The UcS value thus obtained includes a screening effect 
due to the relaxation of the valence-electron density around 
the specific atom. 

The total energies can be basically calculated in density- 
functional formalism 1 15], but an important point is how to 
impose the constraint mentioned above. To this end, we first 
define the d occupation number of the atom / as 
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where ^ 



is a Bloch orbital of a band a with a wave vector 



k, f^j^ is its occupation number, {w/^} are d-type maximally 
localized Wannier functions centered at the atom /, and the 
index /.i specifies the five types of the d orbital. With the def- 
inition of the Njd above, we may write a constrained density 
functional as 



E{Nid) = mm 
p 



|i^[p(r)] - A 
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Here, F[p] is a usual density functional with a total charge 
density p{r) = J2ak l'?^ctfc(''')l^' ^'^'^ A is a Lagrange multi- 
plier. A functional derivative of E{Njd) with respect to the 
Bloch orbital (j>^f^ leads to the following constrained Kohn- 
Sham (KS) equation. 



Hks + -^E \'^if^)i'^if^\ 



^ak - ^ak^ak^ 



(4) 



where Hks is a one-body KS Hamiltonian, and the second 
term in the brackets is an additional potential due to the con- 
straint. 

The calculation for the effective onsite Coulomb interac- 
tion UcS proceeds as follows: We first perform a total-energy 
density-functional calculation with no constraint (A ~ 0) to 
obtain the ground-state energy E{n). Then, we transform the 
resulting Bloch functions {4>^^} into the maximally local- 
ized WF's {ifj^"} used as the input of the subsequent con- 
strained calculations. We then solve the constrained KS equa- 
tion [Eq. (|4}] with a non-zero A, and calculate the occupation 
number Njd in Eq. (|2j using the resulting Bloch functions 
{4>^^} and the WF's above {u'^^y'^}. The A parameter is ad- 
justed so that Nid should be a desired number; i.e., n — 1 
ox n + 1. Finally, we calculate Ucs of Eq. ([3 with the total 
energies thus obtained, E{n), E{n — 1), and E{n + 1). 

We implemented the scheme presented here in Tokyo Ab 
initio Program Package |16] developed by the condensed- 
matter- theory group in the University of Tokyo. With this pro- 
gram, band calculations were performed within the general- 
ized gradient approximation 11711 to density-functional theory, 
using a plane-wave basis set and the Troullier-Martins norm- 
conserving pseudopotentials |18| in the Kleinman-Bylanger 



representation fl^ . The energy cutoff was set to 64 Ry, and 
a 9x9x9 fc-point sampling was employed to represent elec- 
tronic structures of transition metals |20]. A spin polarization 
effect was neglected in the present calculations, and this ef- 
fect seems to be negligible for the evaluation of U^s 1 1 1 1 ■ The 
d-type WF's were constructed with the algorithm of Souza, 
Marzari, and Vanderbilt 1 1 3 1 . To show the validity of repre- 
senting d bands by the resulting WF's, we compare in Fig. 1 
original bands of Fe (solid line) with interpolated bands (dot- 
ted line) obtained by diagonalizing fc-space Hamiltonian ma- 
trix represented by the five d-type WF's. We see that the cal- 
culated interpolated bands reproduce reasonably the original 
narrow d bands, except when sp bands cross the d bands. Con- 
strained calculations were performed for a suppercell contain- 
ing 8 transition-metal atoms to consider a full relaxation of a 
screening charge density to a local charge perturbation. It was 
found that the size of the suppercell is sufficient for obtaining 
a converged f/cff value. 
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FIG. 1: Solid line: Calculated band structure of Fe. Dotted line: 
Interpolated bands obtained from the five d-type Wannier functions. 
Energy window fisll used to compute the Wannier functions was set 
to [— 9.0,+8.5] eV. The zero of the energy scale is at the Fermi en- 
ergy. 

We show in Fig. 2 the calculated effective onsite Coulomb 
interaction U^s (solid line) of the 3c? metals. An importance of 
a screening effect is clear from a comparison between the UcS 
value and a bare Coulomb integral averaged over five d-type 
WF's, 
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f/barc = 5 E 



|w/p(r)|2|u;/^(r')P 



■drdr' . 
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As expected, an unscreened value [/bare (dotted line) is rather 
large relative to UcS- The screening effect makes a signifi- 
cant reduction of the value, thus leading the theoretical values 
to the same energy order as experimental values (dashed | ^ 
and dotted-dashed |6| lines) deduced from a combined use 
of Auger and x-ray photoemission spectroscopy. A more de- 
tailed comparison with the t/cxpt's requires to calculate ab ini- 
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FIG. 2: UcS (solid line), ?7baic (clotted line), and two experimental 
onsite Coulomb interactions (dashed and dotted-dashed lines). The 
experimental values are taken from Ref. |5] for L'^^'-Jpt ™d IS 
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tio Auger and photoemission spectra, but this is not the scope 
of the present study. 

An interesting behavior seen in the figure is that Ucs in- 
creases gradually as an atomic number increases, and this 
trend agrees with the experiments. To our knowledge, there 
is no explanation for the chemical trend observed in [/cfr, so 
we focus on understanding this behavior. If the screening ef- 
fect is neglected, Ucs is Uha,ic, for which it is apparent from 
Eq. that the value of f/barc reflects a feature of the WF it- 
self. More specifically, the observed increasing behavior of 
C/baro results from the trend that the WF shrinks as an early 
metal goes to a late one, which will be also stated in terms 
of the Wannier spread; namely, [/bare is inversely scaled by a 
locahzation length describing a spatial extent of WF 



Q. = 



(6) 



where r is electron coordinates and (r)^ = J r\w i ^{r)]"^ dr , 
etc. 

In contrast to C/barc, the analysis for UcS itself is rather in- 
tractable, because, in that case, in addition to the shrinking 
effect above, a difference in the screening in the metals also 
contributes to the trend in Ucs- In the transition metal, the 
screening results mainly from itinerant sp electrons behaving 
as a free electron, so the effective Coulomb interaction be- 
tween two d electrons is to a first approximation written as 
an interaction between test charges placed in a homogeneous 
electron gas, i.e., as the Yukawa-type interaction e^'<^l'^'^) jr 
with CTo being the Thomas-Fermi screening length ll2lll . From 
this viewpoint, we give a heuristic formula for C/cfr, which is 
assumed to be a one-center integral of the Yukawa potential 



C/eff = 



v'\la) 



r — r' 



x|w/^(r')prfrrfr', 



(7) 



where a in the integral is an adjustable parameter to specify 
a characteristic length of the screening in the system. As the 
value of a becomes larger (i.e., the screening is weakened), 
IJaS approaches to J/barc (unscreened value), so this parame- 
ter is regarded as a measure of the screening strength. We de- 
termined the a parameter so that the integral value [/off should 
coincide with the C/cg obtained from ah initio calculations. By 
plotting a for all the 3d species, we can discuss the trend in 
the screening strength of the transition metals. 

We display in Fig. 3 the resulting dependence of a on 
atomic species (dashed line), together with that of (dotted 
line). The figure clearly shows that the increasing behavior 
of C/cff (solid line) results from the orbital shrink (i.e., the de- 
crease of ri). In contrast, the effective screening length cr is 
almost the same for the listed metals, although there is a small 
jump of cr between Ni and Cu, which reflects the fact that Ni 
has two screening channels of sp and d electrons, while Cu 
has no latter channel because the Cu d bands are fully occu- 
pied. We note that the calculated a agrees reasonably well 
with the Thomas-Fermi screening length ctq (dotted-dashed 
line) of a free electron gas with the same density as the sp- 
electron density of the system i23ll . which indicates a validity 
of our cr estimation based on Eq. 0. The small difference 
between cr and ctq 0.1 A) reflects the presence or absence 
of a screening contribution from intersite d electron transfers. 
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FIG. 3: t/cff (solid line), Wannier spread SI (dotted line), effec- 
tive screening length a (dashed line), and Thomas-Fermi screening 
length (Jo (dotted-dashed line). 

In summary, we have developed a constrained local den- 
sity functional approach with maximally localized Wannier 
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functions, and applied it to the systematic study of the ef- 
fective onsite Coulomb interaction of 3d transition metals. 
This scheme is of practical significance because it is free from 
a limitation on a choice of basis functions, and thus allows 
us to perform the constrained calculations with plane-wave- 
based electronic-structure codes. This property will be partic- 
ularly helpful in a situation where one's interest is in molec- 
ular solids such as BEDT-TTF |24] and solid Ceo iil E3l . 
for which the plane-wave basis is suitable for calculating their 
electronic structures. To see whether the present scheme can 
be indeed exploited in constructing ab initio model Hamil- 
tonians and/or in discussing experimental J7's for the above 
mentioned systems, however, needs future studies. 
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